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The structural, electronic and transport properties of metallocene molecules (MCp2) and isolated or nanotube- 
encapsulated metallocene chains are studied by using a combination of density functional theory and non- 
equilibrium Green's functions. The analysis first discusses the whole series of isolated MCp2 molecules, where 
M = V, Cr, Mn, Fe, Co, Ni, Ru, and Os. The series presents a rich range of electronic and magnetic behaviors due 
to the interplay between the crystal field interaction and Hund's rules, as the occupation of the d shell increases. 
The article then shows how many of these interesting properties can also be seen when MCp2 molecules are 
linked together to form periodic chains. Interestingly, a large portion of these chains display metallic, and even- 
tually magnetic, behavior. These properties may render these systems as useful tools for spintronics applications 
but this is hindered by the lack of mechanical stability of the chains. It is finally argued that encapsulation of the 
chains inside carbon nanotubes, that is exothermic for radii larger than 4.5 A, provides the missing mechanical 
stability and electrical isolation. The structural stability, charge transfer, magnetic and electronic behavior of 
the ensuing chains, as well as the modification of the electrostatic potential in the nanotube wall produced by 
the metallocenes are thoroughly discussed. We argue that the full devices can be characterized by two doped, 
strongly correlated Hubbard models whose mutual hybridization is almost negligible. The charge transferred 
from the chains produces a strong modification of the electrostatic potential in the nanotube walls, which is am- 
plified in case of semiconducting and endothermic nanotubes. The transport properties of isolated metallocenes 
between semi-infinite nanotubes are also analyzed and shown to lead to important changes in the transmission 
coefficients of clean nanotubes for high energies. 

PACS numbers: 85.35.Kt,85.65.+h,71.15.Mb,71.10.Pm 



I. INTRODUCTION 

Carbon nanotubes 1 have been for long proposed as possi- 
ble elements for electronic circuits on the nanoscale. Apart 
from their mechanical stability and small diameter, these de- 
vices can be produced in large quantities and be easily assem- 
bled. Furthermore, they can be shaped in very different geo- 
metric configurations which make them specially suitable to 
join different parts of a nanocircuit and achieve novel switch- 
ing features 2 . On the other hand the electronic properties can 
range from insulating to metallic 3 depending on the radius 
and the chirality. In the absence of defects, electron transport 
falls in the ballistic regime and therefore large carrier mobil- 
ity, high yield and low power consumption is expected. Ad- 
ditional functionality can be achieved by promoting defects 
in the nanotube walls 6 -, or by doping them, n- and /?-doped 
nanotubes can be fabricated by pulling charge in or out of 
the walla 7 ! 8 ! 9 ! 10 ] 11 ! 12 ! 13 ] 14 ] 15 ! 16 . These functionalities enable the 
nanotubes to mimic different parts of a transistof 4 ^ 6 - and even 
to produce devices with negative differential resistance 8 . 

Several methods have been used to produce n- and p-doped 
SWNTs but most of them have the problem that the result- 
ing SWNTs are not stable upon air exposition due to oxida- 
tion processes which tend to /?-dope the carbon wall. A bet- 
ter procedure consists of encapsulating 17 metallic atomsi£, or 
even organic molecules like fullerenes 19 inside the nanotube 
walls. Metallic elements and compounds tend to form contin- 
uous conducting nanowires 17 , while fullerenes produce local 
distortions in the electronic band structure of the nanotube 19 . 
This way, fullerenes can be used to section the nanotube into 
different segments^ -. This local change of the gap should ef- 
fectively produce a series of quantum dots that could local- 



ize electrons 2 -. Trapping single electrons may be a used to 
produce entangled electron pairs and also to tailor interac- 
tions between static and flying qubits 22 . The type of charge 
transfer can be traced back to electron affinity and the ioniza- 
tion potential of the encapsulated molecules. It is recognized 
now that molecules with high electron affinity tend to remove 
charge from the nanotube wall and produce ;?-doped SWNTs, 
whereas molecules with low ionization potential tend to trans- 
fer charge to the wall and produce n-doped SWNT s 23 ' 24 . 

Nanotubes can also be rendered magnetically active 
upon encapsulation of magnetic elements or organometallic 
molecules. These new devices offer the possibility to sepa- 
rate the magnetic material from the environment and hence 
reduce its oxidation. Additionally, encapsulation may give 
rise to new magnetic properties not present in conventional 
magnetic nanowires 25 and even change the mechanical sta- 
bility of the nanotubes 26 . These include, for example, hys- 
teresis loop shifts due to exchange bias 27,28 , changes in the 
complex permittivity and permeability 29 , increased coercitiv- 
ity and Barkhausen jumps 30 , magnetic anisotropy 31 , size and 
number-dependent interactions between catalyst particles 32 
and many other magnetic phenomena. Other fascinating prop- 
erties include the disappearance of magnetism in small clus- 
ters of Ni atoms 34 or the antiferromagnetic behavior shown by 
irorOl 36 . 

The subject-matter of this article is a theoretical proposal 
for the encapsulation of specific metallorganic molecules in 
nanotubes, that combine the two functionalities sketched in 
the former two paragraphs, e. g.: tailored doping and addi- 
tion of magnetic behavior to carbon nanotubes. The proposed 
metallorganic molecules are called metallocene s 38 ' 39 and are 
composed of a metallic transition metal (TM) atom and two 
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cyclopentadyenil (Cp) rings made out of five carbon and five 
hydrogen atoms each. The configuration of both rings is spe- 
cially stable if they loose an electron, thus transforming in 
aromatic rings where the 7r-electrons are delocalized around 
the ring. This produces a negative compound that can bind 
to another Cp ring through a TM atom with a valence of +2 
or higher. The case with a valence of +4, which is specially 
relevant for the fabrication of polymers is not dealt with in the 
present paper. The article instead focuses only on the sim- 
plest metallocene molecule structures, which have only two 
Cp rings bridged by the TM atom. We note that some of these 
molecules have already been suggested as the functional core 
of ingenious electronics devices. For instance, FeCp2 (fer- 
rocene) has recently been proposed as a molecular bridge that 
can rival the conductance of conjugated molecules 40,41 and as 
a molecular switch 43 . A spintronics switch made out of two 
CoCp2 molecules (cobaltocenes), sandwiched by two mag- 
netic leads, has also been proposed 42 . Explicit encapsula- 
tion of cobaltocene in carbon nanotubes was recently real- 
ized by Li and coworkers 37 , who showed that such compos- 
ite devices are stable for nanotube diameters of 9.4 A. For 
thicker nanotubes only a modified version of the cobaltocene 
(bis(ethylcyclopentadienyl)) was observed to remain inside 
the tube. This was possibly due to the fact that cobaltocenes, 
which are smaller than bis(ethylcyclopentadienyl), bonded too 
weakly to the wall of the thicker tubes and were expelled 
from them at the slightest thermal or mechanical fluctuations. 
Li and coworkers also reported a charge transfer from the 
molecule to the nanotube wall and a red shift in the photolu- 
minescence spectra due to the formation of localized impurity 
states, which could make these systems specially suitable to 
produce n-doped nanotubes and modulate the electronic prop- 
erties of the nanotube wall. 

Noting that TMCp2 molecules could be synthesized with 
TM= V, Cr, Mn, Co and Ni, as well as Ru and Os, and that each 
of them display different electronic and magnetic properties, 
we proposed 44 to use them as off-the-shell elements to fab- 
ricate one-dimensional chains of metallocene molecules, en- 
capsulated inside the carbon nanotubes walls. In other words, 
we proposed to fabricate nanoscale wires whose tailored elec- 
tronic, magnetic and conducting properties would be due to 
the transition metal atoms of the metallocene molecules, and 
where the major role of the nanotube would be to actually 
provide mechanical and electric isolation to the metallocene 
core. We envisioned two possible categories of devices. In 
the first class, the nanotube walls would also be metallic, and 
therefore the properties of the device would come from a non- 
conventional interplay between two non-trivial electron sys- 
tems. In the second category, the nanotube walls remain semi- 
conducting, and the analog of a electrical appliance wire on 
the nanoscale would be realized. 

In this article we take further our previous calculations and 
show a host of additional features that can be present in these 
systems. We analyze the structural and magnetic properties of 
isolated metallocene molecules and isolated or encapsulated 
chains. We confirm that the devices display a rich variety 
of electronic properties that can indeed be traced back to the 
diverse filling of the metallocene molecular levels. We find 



that the metallocene chains are very accurate realizations of 
the one-dimensional Hubbard model, with hopping, Hubbard 
U and doping parameters that can be extracted unequivocally 
from the simulations. These chains transfer charge to the nan- 
otube walls, but do not hybridize with them. For the metallic 
nanotube class we therefore find two Hubbard liquids that act 
as reservoirs of each other. On the contrary, for the semi- 
conducting nanotube class, and whenever the encapsulation 
is endothermic (e.g.: when the diameter of the nanotube is 
small), we find that the metallocene molecules produce sig- 
nificant energy barriers on the nanotube walls that could lo- 
calize electrons. We believe that these devices may be useful 
for quantum computation, due to the possibility to realize and 
tailor interactions of localized and flying qubits 22 . 

The outline of the article is as follows: we describe in Sec- 
tion [n] our assembly of theoretical methods. The properties 
of isolated metallocene molecules and chains of metallocenes 
are shown in Section [HI] In Sections [TV] and [V] we focus on 
the properties of chains of metallocenes inside nanotubes and 
the mapping to strongly correlated models. In Section [VT| we 
study the transport properties of isolated metallocenes inside 
carbon nanotubes. We close the article with a section that 
presents and discusses our conclusions. 



II. THEORETICAL METHOD 

The calculations were based on Density Functional 
Theory 45 as implemented in the first-principles code 
SIESTA 46 , which uses norm-conserving pseudopotentials 4 - 7 - to 
get rid of the core electrons and linear combinations of local- 
ized atomic orbitals <f>^(f — d^) to span the valence states 48 . 
The pseudopotentials, which were optimized so that the local 
part was smooth 46 , also included non-linear core corrections 49 
for the transition metal atoms of the MCp2 molecules, with 
core correction radii of 0.70 a. u. These were obtained by 
following the recipe used to generate the pseudopotential in 
iron explained in Ref. |5CJ. The basis set had a double— £ char- 
acter for carbon and hydrogen atoms and double— ( polarized 
for the transition metal atoms. The cutoff radii were deter- 
mined by minimizing the energy for a given reference system, 
i. e.: the H2 molecule for hydrogen, a diamond crystal for car- 
bon, and the bulk phase for the corresponding transition metal 
atom. 

Due to the localized character of the atomic orbitals spe- 
cial care had to be taken to minimize the impact on our results 
of both the basis set superposition error (BSSE) 52 and the ba- 
sis set incompletion error (BSIE) 53,54 . The first error comes 
from the fact that the Hilbert space of the composite system 
is different from the Hilbert space of the each part taken sep- 
arately and therefore the binding energy is smaller due to an 
additional contribution which comes from the the presence of 
more localized states. In order to correct BSSE, it is neces- 
sary then to subtract the additional contribution by calculating 
in separate simulations the energy of each individual system 
in the geometry of the composite system, but inserting ghost 
states at the location of the missing atoms. Then, the binding 
energy of the device is estimated as 
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AE = AE - (E A+GB - E A ) - (E B+GA - Eg) (1) 

where A and B represent the nanotube and the molecule in 
the relaxed configuration of the composite system, GA and 
GB are the corresponding ghost states and AE Q = E A b — 
E Alcl — EB rcl is the binding energy calculated without ghost 
states and relaxing all structures. 

The contribution from the BSIE appears when there is 
large interstitial separation between different parts of a sys- 
tem which is described with localized basis sets. Although 
it is not possible to estimate this contribution it is known to 
have an opposite sign to the BSSE. This implies that the calcu- 
lated values of the BSSE-corrected binding energy are upper 
bounds. We note anyway that we were careful enough to span 
all the interstitial space between the molecule and the nan- 
otube wall so that the states of the molecule and the nanotube 
wall had enough overlap. 

The exchange and correlation energy and potential were 
calculated with the generalized gradient approximation 
(GGA) as parameterized by Perdew, Burke and Ernzerhof 55 , 
since this approximation works much better for transition 
metal elements 50,5 ! . The real space grid used to compute the 
integrals leading to density and Hamiltonian matrix elements 
was calculated by setting a plane wave cutoff of 200 Ry. Fi- 
nally, we took a number of A; -points along the direction of the 
nanotube (z-axis), ranging from 20 to 100, depending on the 
size of the system and on wether the nanotube was semicon- 
ducting or metallic. 

The transport properties were obtained with the ab-initio 
code SMEAGOL 57 , which uses the Hamiltonian provided by 
SIESTA to calculate the retarded Green's function of the scat- 
tering region (M), which in the basis of atomic orbitals is ex- 
pressed as 

Gym^E) = [E+S -H- - Z*(E)}^ M ^ (2) 

where S is the overlap matrix, H is the Hamiltonian, E + = 

E + iS and 

E*(£) = [E+S - H] ML Gf; h (E)[E+S - H] LM (3) 
Sg(£?) = [E+S - H] MR G^ R (E)[E+S - H] RM (4) 

are the self-energies that take into account the coupling to the 
left and right semi-infinite electrodes. Gll(rr) are the surface 
Green's functions. From the retarded Green's function and 
its Hermitian conjugate, the advanced Green's function G A , 
it is possible to obtain the density matrix using the Keldysh 
formalism 58 , so that 



= / dEG< M ^(E) (5) 

where G < is the retarded Green's function, which in general 
is calculated as 



G< M (E) = G(E)l M [T L f L + T R f R ](E)G(E)^ M (6) 

where r L ( R ) = ipwjy — ^ A (R)] are ^ e g amma matrices that 
take into account the strength of the coupling to the leads and 
fh(R)(E) = f(E — //l(r)) are the Fermi distribution func- 
tions of the leads. Such a procedure is necessary to obtain the 
electronic density, n(r) = Yj^v P^A^r - d^ v {f - d v ), 
because the calculation deals with semi-infinite leads and also 
in a general non-equilibrium situation the distribution func- 
tion of the scattering region is not well defined and therefore 
the common diagonalization procedure can not be used. 

Once the self-consistency is achieved the conductances for 
spin-up and spin-down electrons are calculated by using the 
familiar expression^ 

g„{E) = ^tr[r L G A M r R G* M r R M£) (7) 

where a — {f, J.} is the spin index. The current is simply the 
integral of the conductance multiplied by the difference of the 
Fermi distribution functions of the leads. 



in. ISOLATED TMCP 2 MOLECULES AND CHAINS 

As stated above, metallocenes where the transition metal 
belong to the 3d row or the iron column (ruthenocene and 
osmocene) are composed of a metallic atom sandwiched by 
two Cp rings. Such rings can be located just one on top of 
the other, so that the axis that passes through the middle of 
the rings and the metallic atom is G5 and therefore the sym- 
metry is an eclipsed D$h, or rotated 36 degrees relative to 
each other, so that the symmetry is D$d, i- e. staggered, and 
each carbon or hydrogen atom sits on top of one of the edges 
of the other ring. In many cases the same metalloecene can 
have both structures depending on the phase. For instance, 
vanadocene has D^d symmetry in the solid state, but D§h in 
the gas phase. However, we always found in our calculations 
the first configuration to be more stable than the second one 
by an amount that could be as large as 0.03 eV (for the cobal- 
tocene). This difference can be due to the BSIE, that artifi- 
cially decreases the interaction between the rings. We have 
anyhow tested for several of the calculations that both geo- 
metric arrangements of the Cp rings produce almost identical 
results. We are therefore confident that this difference is not 
important for our study. 

The theoretical average distance between carbon atoms in 
the Cp2 rings is 1.46±0.05 A, although it can change slightly 
depending on the type of metallocene. These results compare 
rather well with experiments, which give an average distance 
of 1.44 A for the ferrocene 39 , for example. The distance be- 
tween carbon and hydrogen atoms is equal to 1.10 A. The 
main structural parameter, which is the distance between TM 
and carbon atoms is listed in table UJ along with the distance 
between the TM atom and the plane of the Cp ring. In al- 
most all cases the first type of distance is the same for each of 
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FIG. 1: Densities of States g t of (a) VCp2, (b) FeCp2 and (c) 
CoCp2, projected onto the d-shell orbitals of the 3d-row metal atom, 
i — d x 2_ y 2,d X z,d 3z 2_ r 2,dyz,d xy . Continuous and dashed lines 
represent spin-up and spin-down electrons, respectively. 



the ten carbon atoms in each metallocene. However, the axis 
of the two Cp rings in chromocene and, to a lesser extent, in 
cobaltocene, were found to be tilted, in agreement with exper- 
iments, possibly due to the lifting of degeneracy by a Jahn- 
Teller distortion 59 . The distances in these cases vary between 
2.09 and 2.26 A for CrCp 2 , and between 2.07 and 2. 15 A for 
CoCp2- Table J] shows in these two case the average value. 
As can be seen the agreement between theory and experiment 
is excellent, with the only exception of manganocene. The 
reason of such discrepancy could be due to the fact that the 
distances for MnCp2 correspond to the high spin state while 
we only considered the low spin configuration. 

The differences in the electronic structure of the metal- 
locenes can be traced back to the progressive filling of the 
d-shell of the TM atom, that we also show in table U These 
electronic structures are determined by the interplay among 
electron filling, Hund's rules and the crystal field interaction. 
To illustrate this, we show in Fig. (Q~|i the density of states 
projected onto the five d-orbitals of the TM atom, of VCp2, 
FeCp2 and CoCp2- We find, first, that the d xz and d yz or- 
bitals feel a different crystal field than the rest of orbitals, 
since they are oriented towards the carbon atoms in the Cp 
rings. They also overlap strongly with the carbon orbitals 
and, consequently, their energy levels are split into bonding 
and antibonding states. These two orbitals are degenerate in 
energy for all metallocenes, see Fig. (H](a), (b)). Chromocene 
and cobaltocene are an exception since, due to the Jahn-Teller 
distortion the axis of the molecule is bent. Consequently, the 
degeneracy of the d xz and d yz orbitals is lifted, as shown in 
Fig. (Q](c)). The energy levels of the d x 2_ y 2, d xy and d^ r 2_ z 2 
are almost degenerate and lie lower in energy because they are 
not oriented towards the carbon atoms. Additionally, d x 2_ y 2 
and d xy orbitals overlap a little with carbon orbitals, while the 
d 3r 2_ z 2 does not. Iron d-orbitals in FeCp2 thereby arrange in 



TABLE I: Theoretical and experimental distances between the TM 
and the carbon atoms, between the carbon atoms and the Cp2 rings, 
and number of unpaired electrons. Experimental values from Ref. 
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a closed shell structure; consequently FeCp 2 has no unpaired 
electrons and is non-magnetic. The remaining metallocenes 
display a sort of mirror symmetry about the closed-shell fer- 
rocene. Notice that the exchange splittings, shown in Fig. (Q}, 
are proportional to the number of unpaired electrons. It is 
also interesting to note that the geometric arrangement is such 
that the orbital angular moment is quenched, and therefore the 
measured values of the magnetization almost coincide with 
the spin moments. 

The distribution of electrons in the d orbitals can also ex- 
plain why the distance between the Cp ring and TM atoms 
has a minimum for ferrocene. We note that the same be- 
havior is found for the lattice constant of periodic isolated 
metallocene chains 44 , where the C5 axis of the molecule and 
of the chain lie parallel to each other. The electronic band 
structure of isolated VCp2 and CoCp2 chains is plotted in 
Fig. (O. Notice that VCp2 chains are semiconducting with 
a large band gap. Moreover, their bands are essentially flat. 
This is expected, since metallocene molecules are fairly stable 
and should show small overlaps among each other. The band 
structure of CoCp2 chains is more interesting. First, the band 
width is larger -of about a quarter of an electron volt- and, 
second CoCp2 chains are semiconducting, as expected, but 
the gap is fairly small. Both chains are strong ferromagnets, 
since their exchange splittings are such that up- and down-spin 
bands do not overlap in energy. CoCp2 chains are therefore 
more promising for fabrication of transport devices, and this 
paper shall mostly focus on them henceforth. 

We remind that CoCp2 molecules possess a bent axis, that 
leads to a lifting of degeneracy of the molecular levels asso- 
ciated to the d xz and d yz orbitals. Fig. (|2](b)) shows that the 
molecular crystal field also lifts the degeneracy of the one- 
dimensional bands closest to the Fermi energy, that are asso- 
ciated with these two orbitals. We note that the splitting is of 
about 0.5 eV. The figure also shows that these two bands are 
exchange-splitted by an amount of 1.3 eV for the d xz and 0.4 
eV for the d yz . We take now the spin-up and down d xz bands 
and realize that they may be interpreted as the mean-field fer- 
romagnetic solution of a one-dimensional nearest-neighbors 
Hubbard model at half-filling, with a hopping integral of 0.06 
eV (bandwidth of 0.25 eV), provided that there are no inter- 
band excitations (e.g.: electrons from the filled majority-spin, 
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FIG. 2: Band structure around the Fermi level for chains of 
vanadocene ((al) and (a2) for majority and minority spin, respec- 
tively) and cobaltocene ((bl) and (b2)) calculated at a lattice con- 
stant of 7.5 A. Circles, squares, diamonds, triangles up and triangles 
down represent d yz , d xz (e* 9 orbitals), d x 2_ y 2, d xy (ei g orbitals) 
and d 3z 2_ r 2 (a' lg orbital), respectively. 



d xz band can not be excited to neither of the unfilled majority 
or minority spin, d yz bands). The coulomb integral U can then 
be extracted by noting that the spin-splitting is equal to Un, 
where n = 1 is the filling fraction. Hence U is equal to 1.3 
eV. But the ground state of a half-filled Hubbard model is not 
a ferromagnetic metal, but an insulating antiferromagnet. We 
have therefore simulated helical arrangements of spin s 50 ' 51 . 
We show in Fig. ([3] (a)) the energy of these spin spirals as 
a function of the wave-vector q of the helix. The figure shows 
that the ground state is indeed an antiferromagnet. This fact 
further confirms our belief that these chains are perfect play- 
grounds of Hubbard physics. A final word of caution should 
however be added, since according to the former arguments, 
the empty spin-up and -down d yz bands should not be splitted 
at all. 



IV. CHAINS OF TMCP 2 MOLECULES INSIDE CARBON 
NANOTUBES 

A. Motivation and notation 

Isolated metallocene chains clearly lack mechanical stabil- 
ity. We therefore propose that encapsulating them inside car- 
bon nanotubes should provide the wanted isolation and sta- 
bilization. Inclusion of organic molecules inside carbon nan- 
otubes entails however a series of stringent requirements. En- 
capsulation of individual molecules can surely be achieved if 
the chemical process is exothermic. But assumed this con- 
dition is fulfilled, the energy barriers for movement of the 
molecules along the axis of the nanotube should be tuned to 
be not too high to completely prevent this motion, nor too low 
to allow the molecules to displace freely. Notice that in this 
last case molecules should leave the tube under slight mechan- 
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FIG. 3: Energy as a function of the spiral pitch vector for (a) an iso- 
lated chain of cobaltocenes and (b) a chain of cobaltocenes inside a 
(7,7) SWNT grown along the C5 symmetry axis, q = corresponds 
to the ferromagnetic configuration and q — 1 (in units of n/a) to the 
antiferromagnetic configuration. 

ical shake-ups. Chain formation from the individual encapsu- 
lated molecules is helped by the own binding energy of the 
chain. This additional source of stabilization energy should 
link molecules together, and counteract the possible negative 
effect of the detuning of the energy barriers. Endothermic en- 
capsulation of fullerenes has also been achieved^ .. We there- 
fore believe that individual metallocene molecules may be en- 
capsulated by endothermic processes. 

Encapsulation produces mechanical and electrical isolation 
of the metallocene chain. It also dopes the nanotube and pro- 
duces local modifications of its band structure 20 . These local 
modifications can be seen as electrostatic wells that trap elec- 
trons. We expect these traps to be deepest for endothermic en- 
capsulation since in such a case the geometry would enhance 
the interaction between molecules and nanotube walls. Notice 
also that the length of the traps along the nanotube should be 
small for metallic nanotubes since in this case the screening 
length is pretty short. On the contrary, the length of the traps 
should be much larger for semiconducting nanotubes. 

We wish to stress the importance of commensurating the 
lattice constant of the isolated chain (7.50 A) with that of 
the nanotube, so that molecules accommodate well inside the 
tube and the chain is neither compressed nor strained. If the 
metallocene molecules are too close, electron clouds of two 
neighboring molecules will overlap too much and the full 
device is rendered endothermic. If they are too separated, 
then their wave functions won't overlap appreciably and the 
molecules can be considered as isolated. Since the unit cell of 
a armchair nanotube has a lattice constant of 2.50 A, a met- 
allocene may be well accomodated every three nanotube unit 
cells. The lattice constant of zigzag naotubes is 4.32 A, so 
placing a molecule every two nanotube unit cells renders a 
stretched chain. Chiral tubes like the (12,4) or (15,5) have 
a lattice constant of 15.59 A, and therefore one can safely 
place two metallocenes inside each of them to form a some- 
waht stretched chain. To avoid confusion, we use the follow- 
ing notation to denote the basic unit cell of the simulations: 
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FIG. 4: Cohesive energy for the cases N = 3 (circles) and N = 4 
(squares) in the parallel (continuous line) and perpendicular (broken 
line) configurations, calculated without (a) and with (b) BSSE cor- 
rection. 




FIG. 5: Band structure of (a) an isolated cobaltocene chain (lattice 
constant = 7.5 A), (b) a cobaltocene inside 3 unit cells of a (7,7) 
nanotube (the same separation between molecules as before) and (c) 
a (7,7) nanotube, for spin-up and spin-down electrons. 



#TMCp2@-/V(n, n'), where # is the number of metallocenes 
(TMCp 2 ) that are included in N unit cells of a (n, n') SWNT. 
We have also placed the molecules in two possible orienta- 
tions, which we call parallel or perpendicular depending on 
whether the C5 axis of the nanotube is aligned parallel or per- 
pendicular to the nanotube axis. 



B. Metallic nanotubes 

We have found that the BSSE artificially enhances the sta- 
bility of the encapsulated chains, so that BSSE corrections are 
essential to get right answers. This is explicitly shown in Fig. 
dU for CoCp2 chains, encapsulated in armchair tubes of dif- 
ferent diameter. Notice that the most stable configuration cor- 
responds to a (8,8), or (7,7) nanotubes (radii of 5.4 or 4.7 A, 
respectively) depending on whether BSSE is corrected or not. 
Hence, we shall only discuss our results with BSSE correc- 
tion, Fig. (0](b)). We note that encapsulation is exothermic 
for armchair (n, n) tubes, with n equal or larger than 7 (radii 
equal or larger than 4,7 A). The figure predicts that the highest 
rate of encapsulation should be achieved for (8,8) tubes, with 
smaller rates for (7,7), (9,9), (10,10) and so on. But the TEM 
images of Li et al? 1 only show encapsulation for radii simi- 
lar to that of the (7,7) or slightly larger, and therefore predict 
a very selective encapsulation. We believe this is due to the 
small energy barriers that a metallocene would find when it 
moves along the axis of thicker tubes. Since experiments are 
performed at room temperature and samples are transported 
between different places, small barriers wouldn't prevent the 
metallocene to escape when mechanical and thermal pertur- 
bations increase its mobility. 

Before switching to the analysis of the electronic structure 
we shall further analyze the results presented in Fig. (0](b)). 
The figure demonstrates that molecules fit inside (7,7) tubes 
only if their axes are aligned. Misaligned molecules can be 
placed inside thicker tubes, but with a pretty large energy 
cost. We shall therefore only analyze aligned metallocene 



molecules in what follows. According to the same figure met- 
allocenes can also lower significantly their energy by sticking 
close enough, which means that the most stable configuration 
corresponds to a molecular chain inside the nanotube. This 
is also confirmed experimentally by the TEM images of Li et 
al.—, where the cobaltocenes show a clear tendency to appear 
together. 

We turn now to discuss the electronic properties of the pre- 
ferred CoCp2@3(7,7) devices. The first remarkable property 
of these systems is a charge transfer of 0.6 e~ from the cobal- 
tocene to the nanotube, which is due to the high electron affin- 
ity of the latter 24 . Such a charge transfer, which is responsi- 
ble for the reduction of the magnetic moment of the molecule 
from 1 to 0.4 fis, can also be inferred from the shift in the 
position of the spin-up d xz band when the chain is inserted 
in the tube. This shift can be appreciated by comparing the 
band structures of isolated and encapsulated chains, that we 
depict in Figs. d2](a), (b)). This transfer makes the valence 
bands of the nanotube move down, as can be seen in Figs. 
(0 (b), (c)). On the other hand the crystal field splitting of 
the Co d bands decreases from about 0.5 to 0.2 eV. This ef- 
fect is mainly due to the change in the atomic positions of 
the encapsulated cobaltocene. Indeed, by plotting the bands 
of an isolated chain calculated using the coordinates of the 
cobaltocene in the CoCp2@3(7,7) structure we find that the 
change in the crystalline field is very similar. The screen- 
ing and charge transfer also produces a large reduction of the 
Stoner splitting from 1.3 to 0.5 eV. Taking the mean field esti- 
mate for this splitting, AE m nU = 0.4 U, gives a Hubbard 
U of 1 .3 eV again. We therefore infer that screening effects 
produced by the nanotube do not reduce the Coulomb interac- 
tion in these chains. 

The filling of the metallocene chains modifies the mean- 
field ground state, that was insulating and antiferromagnet in 
vacuum, to an helical state with a small pitch vector very close 
to the ferromagnetic solution, as can be appreciated in Fig. 
(0(b)). The small oscillations in the energy are due to the 
application of the spiral rotation to the nanotube states. 



7 



(a) (b) (c) 




r x rr x rr x r 



FIG. 6: Band structure of (a) CoCp 2 @2(1 1,0) (b) 2CoCp 2 @(12,4) 
and (c) 2CoCp2 @(15,5) for spin up and spin down electrons. 



C. Semiconducting nanotubes 

We extend now the same discussion to the case of semi- 
conducting nanotubes. The evolution of the binding energy 
and charge transfer to the nanotube wall as a function of the 
nanotube radius is very similar to the armchair case. Zigzag 
nanotubes start being exothermic for a radius of ~ 4.7, which 
corresponds to the (12,0) tube. The chiral nanotubes (12,4) 
and (15,5), with radii of 5.7 and 7.2 A respectively, have also 
been found to be exothermic. Furthermore, we have found 
that the charge transfer decreases as the radius of the nanotube 
increases, as in armchair tubes. 

We show in Fig. (|6]i the band structures of CoCp2@2(l 1,0), 
2CoCp 2 @(12,4) and 2CoCp 2 @(15,5) as representative of 
stretched and reasonably commensurated chains. The most 
noticeable feature in the first case, which corresponds to an 
endothermic encapsulation, is the severe reduction of the crys- 
talline field between the d yz and d xz bands as a consequence 
of the change in the atomic positions produced by the strong 
interaction with the nanotube wall. Note also that these bands 
are almost non-dispersive due to the large separation between 
metallocenes along the nanotube axis, which are placed about 
8.64 Aapart. The chiral cases are very similar, despite the 
large difference between their radii, although, as expected, the 
charge transfer is larger in the (12,4). In both cases the cobal- 
tocene bands are clearly dispersive and show a crystalline field 
splitting which in the (15,5) is slightly larger due to the smaller 
interaction with the nanotube wall. Note that bands associ- 
ated wit the cobaltocene are folded due to the presence of two 
molecules in the unit cell. 

The Hubbard U can not be cleanly estimated in the en- 
dothermic (1 1,0) encapsulation due the strong modification of 
the cobaltocene bands and the hybridization with the nanobute 
states. However, in the chiral nanotubes by taking into ac- 
count the exchange splittings, which are 0.72 and 0.87 eV in 
the (12,4) and (15,5), respectively, and the occupations of each 
spin up d xz level, which are 0.57 and 0.67 electrons, respec- 
tively, U turns out to be again 1.3 eV, in agreement with the 
armchair nanotube. This demonstrates the robustness of the 
calculation. 



V. MAPPING TO STRONGLY CORRELATED MODELS 
A. Physics at the chain 

The discussion following Figs. (|2](b)) and 0(a), (b)) high- 
lighted the similarities of the band structure, ground state and 
magnetism of isolated or encapsulated metallocene chains to 
strongly correlated Hubbard chains at diverse doping levels. 
We wish to take a closer look to this delicate issue now. 

Indeed, an isolated metallocene chain may be viewed as a 
series of magnetic impurities stacked on top of each other, 
whose wave-functions have little overlap. The result of the 
process of encapsulation may therefore be a sort of Kondo 
chain, where the metallocene chain acts a an impurity lattice 
that hybridizes with the electronic (Luttinger) liquid residing 
in the host nanotube. 

The parameters that tune the model are therefore the hy- 
bridizations among metallocene molecules Vmm, or between 
metallocenes and the nanotube Vmn- We note that Vmm 1S 
rather small for many of the isolated or encapsulated metal- 
locene chains since these molecules are fairly stable and there- 
fore do not like to share electrons. Exceptions are CoCp2 and 
NiCp2, whose bandwidth are sizeable. We shall therefore con- 
fine the discussion below again to cobaltocene chains. In or- 
der to estimate the size of the nanotube-molecule hybridiza- 
tion integral Vmn, we first look into the band structure of the 
composite device in Fig. ((5] (b)). The figure shows that the 
most important effect of encapsulation is a transfer of charge 
that shifts the metallocene bands upwards and the nanotube 
bands downwards, and, associated with it and with a change 
of atomic positions, a reduction of crystal field and exchange 
splittings. Nanotube and metallocene bands close to the Fermi 
level cross each other with no apparent hybridization. This 
fact indicates that Vmn is essentially zero. Additional confir- 
mation comes from the fact that nanotube bands in the encap- 
sulated device are not spin-splitted and therefore the nanotube 
does not magnetize. 

We cross-check this estimate of the hybridization by tak- 
ing a look at the projected density of states (PDOS) of the 
different molecule and nanotube atoms, which is plotted in 
Figs. (0, ([8]) and (0 for the cases of an isolated CoCp2 and a 
(7,7) nanotube, a CoCp 2 @3(7,7) chain and a CoCp 2 @4(7,7) 
chain, respectively. All of these devices have been simulated 
for the same lattice constant, which is that of the nanotube. 
Comparison of Figs. (O and I©/© confirms our view that 
the main effect of encapsulation is a transfer of charge and a 
structural modification that brings about a reduction of crys- 
tal field and exchange splittings. Furthermore, the last panel 
in Figs. ([H) and (O shows no noticeable splitting of the spin- 
resolved DOS projected on the carbon atoms of the nanotube. 

We further notice that the shape of the DOS of cobalt atoms 
does not change too much. The only noticeable changes ap- 
pear in the energy range [-8,-6] eV below the Fermi energy in 
the DOS projected on the hydrogen and carbon atoms. What 
remains to be seen is whether these modification are due to 
structural deformation of the molecules and nanotube walls 
after encapsulation, or to hybridization between them. We 
note that the sharp resonances associated to hydrogen atoms 
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FIG. 7: From top to bottom, projected densities of states on the cobal- 
tocene Co, C and H atoms for a Cobaltocene chain at its equilibrium 
distance, and a clean (7,7) SWNT. Solid and dashed lines indicate 
spin up and down components, respectively. 
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FIG. 8: From top to bottom, projected densities of states on the cobal- 
tocene Co, C and H atoms and the naotube walls in a CoCp2 @ 3(7,7) 
device. Solid and dashed lines indicate spin up and down compo- 
nents, respectively. 



B. Physics at the walls 

We now look at the change in the electrostatic potential at 
the nanotube walls created by the combined effect of charge 
transfer, structural modification, and hybridization with hy- 
drogen atoms. We expect that the potential will show well- 
defined wells around the positions where the metallocene 
molecules are placed, that will be better or worse screened 
depending on the metallice or semiconducting character of 
the walls. These wells, that should trap electron s 20 ! 21 , will 
form a periodic chain with the same lattice constant as the 
underlying metallocene chain. A good estimate of the shape, 
length and depth of these wells, can be provided by calcu- 
lating the electrostatic potentials at the walls of the nanotube 
of a CoCp2@-/V(n, n') device with and without metallocene 
molecules, and substracting one from the other, AV(z) — 
V H ,with{z) - V H , without {z). We note that this potential 
should have axial symmetry and hence only depend on the 
z-coordinate. 

We have checked that AV^ is rather small for metallic de- 
vices since, as expected, the valence electrons at the walls 
screen the charge perturbation efficiently. We have subse- 
quently looked into the semiconducting CoCp2@iV(l 1,0) de- 
vice (radius of 4.4 A), whose band structure is plotted in Fig. 
(0 (a)). In this case the encapsulation is endothermic and 
therefore the modification of the potential is expected to be 
further enhanced as a consequence of the strong interaction 
between the metallocene and the nanotube wall. Indeed we 
have found that the potential barrier is rather large as can be 
seen in Fig. ( fTUl ), where we plot Ay as a function of the axial 
coordinate for a device where N — 2 and for another where 
TV = 4. We note that we have placed one molecule at the ori- 
gin of the z-coordinate. We inmediately note the periodicity 
in AF for both N = 2 (8.6 A) and N = 4 (17.2 A). The 
potential has two symmetric minima, placed at the positions 
of the close-by C5 rings and a maximum in between them, 
that corresponds to the position of the cobalt atom, which is 
farther apart. The absolute maximum in AV(z) appears at 
the middle point between two molecules. The height of the 
well depends on both the lattice constant of the chain and the 
screening length. The figure shows that AV(z) has enough 
room to reach zero when N = 4, allowing us to estimate a 
value of about 6-7 A for the screening length. For N = 2, 
on the contrary, the distance between wells (8.6 A) is not long 
enough to allow the complete healing of the potential. Ac- 
cordingly, the depth of the well decreases from 0.2 to 01 eV. 



of the cobaltocene chain are clearly broadened after encapsu- 
lation. We therefore believe that Hydrogen atoms, and partly 
carbon atoms, do hybridize. This hybridization does have no 
implications in the electronic properties of the devices, since 
they appear well away the Fermi level. They provide a glue 
between molecules and tube, that hinders the axial motion of 
the molecule. 



VI. TRANSPORT PROPERTIES OF ISOLATED 
METALLOCENES INSIDE CARBON NANOTUBES 

The case where cobaltocenes form a chain inside the 
nanotube^ can lead to magnetorresistive ratios of 20% due 
to the suppression of the spin channel that goes through the 
molecular chain. A different situation however is that pro- 
duced when a cobaltocene sits in between two infinite nan- 
otube electrodes, as shown in Fig. ( fTTT ). 

The conductance of such a system is plotted in Fig. ( fl~2l for 
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FIG. 9: From top to bottom, projected densities of states on the cobal- 
tocene Co, C and H atoms and the naotube walls in a CoCp2 @4(7,7) 
device. Solid and dashed lines indicate spin up and down compo- 
nents, respectively. 



three different metallocenes, namely ferrocene, cobaltocene 
and nickelocene in the parallel configuration. These are repre- 
sentative cases since the other metallocenes, with the possible 
exception of the ruthenocene and the osmocene, have smaller 
dimensions and are expected to produce less interaction with 
the nanotube wall. The transmission of the clean nanotube 
is also plotted for comparison. As can be seen the net effect 
produced on the transmission coefficients by isolated metal- 
locenes is not very large and the transmission at the Fermi 
level is not affected. Notice that in the case of the magnetic 
molecules the transmission is exactly the same for spin-up and 
spin-down electron due to the fact that the carbon nanotube 
does not acquire any spin polarization from the charge that is 
transferred from the metallocene. The only differences with 
a clean nanotube appear above and below the fermi level and 
have the effect of rounding the transmission coefficients where 
the number of channels increases. This effect is mainly due to 
the structural deformation induced by the metallocene rather 
than to any electronic influence, although it could also be due 
to the charge transfer. The change seems to be very small for 
the ferrocene impurity whereas the cobaltocene produces the 
largest modifications. 

We therefore predict that there should be no change in the 
transport measurements of isolated metallocenes inside nan- 
otubes with radii larger 4.5 A, unless the applied voltage is 
large enough (> 2 V) to reach the regions with large number 
of channels. In such a case the presence of scattering would 
decrease the measured current and conductance and this effect 
would be clearly seen in the second derivative of the current 
with respect to the voltage. It is also worth mentioning that 
such and effect would be more pronounced in nanotubes with 
small diameter but the encapsulation of metallocenes would 
be difficult to realize due to the endothermic character, which 
would probably produce defects that would have a larger in- 
fluence on the transport properties. 
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FIG. 10: AV as a function of the z-coordinate for CoCp 2 @2(ll,0) 
(continuous line) and CoCp2@4(ll,0) (dashed line) devices. 




FIG. 1 1 : Unit cell used to calculate the transport properties of a clean 
nanotube with a metallocene impurity. The atoms of the metallocene 
have been highlighted for clarity. 



VII. CONCLUSIONS 

We have shown that density functional theory in combina- 
tion with pseudopotentials and linear combinations of atomic 
orbitals provides a very accurate description of the electronic 
and structural properties of metallocenes, whose results are in 
excellent agreement with experiments. The electronic config- 
uration of metallocenes, which is mainly determined by the 
interplay between the crystalline field splitting and Hund's 
rules on the d bands, shows a kind of mirror symmetry with 
respect to the column of iron due to the number of unpaired 
electrons. The filling of the d bands and the shape and sym- 
metry of the orbitals involved can also explain the distribution 
of the structural parameters as a function of the repulsion be- 
tween the orbitals of the cyclopentadienyl rings and those of 
the transition metal atom. As a consequence the metallocenes 
become bigger as the number of electrons of the metallic atom 
increases. 
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FIG. 12: Transmission coefficients of isolated ferrocenes, cobal- 
tocenes and nickelocenes in the parallel configuration between semi- 
infinite (7,7) SWNT leads. For comparison it is also shown the trans- 
mission coefficients of a clean (7,7) SWNT. 

The electronic configuration of metallocene chains does not 
differ too much from that of the isolated metallocenes, the 
main difference being the finite dispersion of the states around 
the Fermi level, which make them conducting. Such states, 
which in the cobaltocene are half-filled are responsible for the 
magnetic configuration of the chains, which we demonstrated 
to be antiferromagnetic by calculating the energy of a series of 
helical configurations with the help of the generalized Bloch's 
theorem. The spin polarization of those states makes these 
systems perfect candidates for spintronics applications. How- 
ever, it is necessary to look for possible methods to stabilize 
them without destroying such interesting properties. 

The most obvious way of producing stable metallocene 
chains is by placing them inside a carbon nanotube, which 
would act as an external coating that would protect them me- 
chanically and electrically. Such a system was found to be 
exothermic for radii larger than 4.5 A, in agreement with re- 
cent experiments. The encapsulated metallocenes produce a 
transfer to the nanotube wall and the carbon nanotube be- 



comes n-doped. This effect is clearly seen when the band 
structure or the density of states of the composite system is 
analyzed. Such charge transfer changes also the most stable 
magnetic configuration which moves close to the ferromag- 
netic point with a small spiral pitch vector, although the dif- 
ference in energy is small enough to point all spins up with 
the application of a magnetic field. By closely analyzing the 
projected density of states it is also possible to study the inter- 
action of the molecular magnetic moment with the conduction 
electrons of the carbon nanotube. 

Our mean-field calculations show that CoCp2 chains be- 
have as strongly correlated electron liquids, that can be per- 
fectly well described by a half-filled Hubbard model. We fur- 
ther note that the hopping t and Coulomb U parameters can 
be univocally extracted from our simulations. We also show 
that encapsulation only produces a charge transfer that carries 
the Hubbard liquid away from half-filling. 

Looking into the effects on the nanotube walls, we have 
shown that encapsulation inside semiconducting nanotube 
produces deep wells in the electrostatic potential that should 
act as quantum dots that trap electrons. 

Finally, we also studied the transport characteristics of met- 
allocenes between semi-infinite clean nanotubes. We found 
that the transmission coefficients do not change too much with 
respect to the transmission coefficients of a clean nanotube, 
the only difference being the rounding of the steps where the 
number of channels increases. This would give rise to changes 
in the second derivative of the current with respect to the volt- 
age for large bias voltages. 
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